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We study the linear response to an external perturbation of a renewal process, in an aging condi- 
tion that, with no perturbation, would yield super-diffusion. We use the phenomenological approach 
to the linear response adopted in earlier work of other groups, and we find that aging may have 
the effect of annihilating any sign of coherent response to harmonic perturbation. We also derive 
the linear response using dynamic arguments and we find a coherent response, although with an 
intensity dying out very slowly. In the case of a step-like perturbation the dynamic arguments 
' yield in the long-time limit a steady signal whose intensity may be significantly smaller than the 

^ | phenomenological approach prediction. 
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I. INTRODUCTION 



The problem of linear response of a dynamic system to an external perturbation is one topic of general interest in 
the field of statistical thermodynamics. Many authors |2| have studied the linear response emergence from microscopic 
chaos, and thus from a condition that is expected to generate ordinary statistical physics Q. In the last few years 
an increasing number of investigators have been addressing the interesting issue of the response to perturbations of 
dynamical systems departing from the condition of ordinary statistical thermodynamics 0,1H . A great attention has 
been devoted to the case of sub-diffusion that makes it possible to fulfill the Einstein relation [j| even if the dynamic 
system generates anomalous diffusion |6(. 
I ' It is important to point out that the same dynamic generator of intermittence yields either super- or sub-diffusion 
. according to whether the long sojourn times refer to the velocity or to the position state, respectively. Let us consider 
t-H ' a regular one-dimensional lattice, with the fixed distance a between the nearest neighbors, and a particle jumping 
from one to another site of this lattice. Let us imagine that the particle sojourns for a long time r in one site and 
that at the end of this sojourn makes a jump from the position na, with n being an integer number ranging from 
-co to +oo, either to the position (n + l)a or to the position (n — l)a, according to the coin-tossing prescription. 
This condition can be realized with a dynamic generator of events, occurring at times to = 0, ti,t^, ■ ■ ■ • • • , with 
the waiting times Ti = — ti corresponding to a given distribution i^{t). Notice that a single jump involves the 
occurrence of two events. One event is the drawing of a number r from the distribution ip{r), and the other event is 
the coin tossing. It is worthwhile to stress that in this article we refer to time t as an integer number. The adoption 
of the continuous time picture is a good approximation made possible by the choice of conditions that involve t > 1. 

With the same dynamic model and the same coin-tossing prescription, we can generate a different random walk 
process, producing super-diffusion. Here the random walker moves with a velocity of fixed modulus W , whose sign 
is determined by two events occurring at the times ti. As in the earlier example, the two events are the drawing of 
the time Ti from the distribution ip( T ) an d the coin tossing that fixes the velocity sign of the walker between time 
U and time ij+i. This corresponds to defining two states, the state |1 >, with velocity W, and the state |2 >, with 
velocity — W. We refer to the time interval Ti = ti + \ — ti as waiting time. Of course, the walker can repeatedly jump 
to the right (left), in the first case, or maintain the same positive (negative) velocity sign, in the second case. The 
two processes are different, insofar as in the first case, the larger the sojourn time the slower the diffusion process, 
while in the second the opposite condition applies: the larger the sojourn times the faster the diffusion process. In 
this paper we adopt the velocity rather than the position picture, the reason for this choice being that, as we shall 
see, it establishes a connection with the interesting phenomenon of non-Poisson stochastic resonance Q. We shall be 
referring throughout to the wide time regions between ti and tj+i, where no unpredictable events occur, as quiescent 
regions. We shall adopt a dynamic model inspired to turbulence and to the Manneville map where the quiescent 
regions are usually denoted as laminar regions. This dynamic model 0,^3 yields the waiting time distribution tp{r) 
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with the inverse power law form: 



^) = (H~l) ¥TT y I , (1) 



with \x > 1 and 

T > 1. (2) 

The condition of Eq. (0) is essential to adopt, whenever convenient, the earlier mentioned continuum time treatment. 
The probability that no event occurs until time t is 

dtW) = [ TT ^) • (3) 

This function is called survival probability. 

When (i > 2, the mean value of these waiting times is given by 

<T> =(^)' (4) 

The condition fi = 2, where this mean value diverges is the border between the region /x > 2, compatible with the 
existence of an infinitely aged condition and the region /i < 2, which is characterized by perennial aging and 
ergodicity breakdown [l2j. In this paper we focus our attention on the condition of perennial aging fi < 2. 

For a detailed discussion of renewal aging we refer the reader to the recent literature on this subject, for instance, 
Refs. [T^. Il4|. The proper quantitative treatment of aging requires preparation, namely, that the time origin t = 
corresponds to an event occurrence. Aging means that the waiting time distribution V-K 7 ") an d the corresponding 
survival probability ^(r) depend on the time at which the observation process begins. We define with the symbol 
tp(t,t') the probability of meeting the first event at time t when the observation begins at time t' < t and the system 
is prepared at time t = 0. ^(t, £') denotes the corresponding survival probability. The functions ^j(t,t') and t') 
refer to a condition where t = fa, namely, a time at which an event occurs, and t' is a time at which no event occurs 
and observation begins. Of course, the Poisson conditions ip{t, t') = ip(t — t') and &(t, t') — ^(t — t') are violated 
|l3l 1 1 ll] but in the specific case where t' = fa—i, namely, when observation begins with an event occurrence. More 
precisely, the functions i^(t) and W(t) refer to the condition r = fa — fa-i, with fa, denoting, according to the earlier 
definitions, times at which events occur. To discuss the response of renewal non-Poisson systems to perturbations we 
study their behavior in the presence of a time dependent perturbation, a condition that forces us to introduce also 
the function tp(t\t'). This is the conditional probability that at time t an event occurs, given the condition that the 
earlier event occurs at time t' . In the limiting case of extremely weak perturbation it is expected that 

mt') = ^-t')- (5) 

In fact, in the unperturbed case, as earlier pointed out, the waiting time distribution is obtained by recording the 
distance between an event and the next, insofar as any time can be selected as a time origin, if it corresponds to 
an event occurrence. Renewal aging does not have anything to do with the system being driven by physical rules 
changing with time, and it is only a consequence of making observation and preparation at different times. This 
condition, and the equality of Eq. (jSJ) with it, is violated by external perturbation. Let us note furthermore that in 
the literature the symbols $ and are usually adopted to denote a probability and a probability density, respectively. 
Due to the adoption of a discrete time representation in this article both symbols denote probabilities, insofar as, for 
instance, ^>(r)dr = ^(t), as a consequence of the fact that for the integration time step we have dr = 1. 

The problem of the response of sub-diffusion to external perturbation in an aging condition has been recently 
addressed by Sokolov, Blumen and Klafter and by Barkai and Cheng 0]. These authors have adopted the 
same phenomenological model as Bertin and Bouchaud [l7j . In this model the external perturbation, which does 
not affect the sojourn time duration, influences the choice between either jumping to the right or to the left, at the 
end of a sojourn: In the unperturbed case there is no bias, this choice being done with a fair coin, and the effect 
of perturbation is to turn this fair coin into an unfair coin, producing a bias at the moment of the choice between 
the left and the right nearest-neighbor site. The phenomenological method will be extended to the velocity picture. 
However, in addition to adopting the velocity picture, and thus the physical condition of super-diffusion rather than 
sub-diffusion, we plan to go beyond the phenomenological model. For this reason, we shall also study the case when 
the perturbation-induced bias is generated through the influence that the external perturbation exerts on the waiting 
time duration. Although this is in principle a hard problem, we find an analytical solution in the limiting case of weak 
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perturbation: moving along the lines of the linear response method yj, we express the response to perturbation in 
terms of the system's unperturbed statistical properties. In the aging regime of this paper the dynamic model yields 
results that are significantly different from those produced by the phenomenological model, and identical to those of 
Ref. |l8j , where the dynamical model was already used although with no demonstration [l9l | . It is important to stress 
that also the adoption of the phenomenological model allows us to obtain original results, insofar as it deals with the 
response of super-diffusion, rather than sub-diffusion [TH llfij. to an external perturbation. 

The outline of this paper is as follows. In Section[n]we illustrate a dynamic model, generating, in the unperturbed 
case, a non-Poisson intermittent signal. We shall show that an external perturbation has the effect of perturbing the 
distribution of waiting times in both velocity states, and we shall discuss how to convert this property into a bias. 
Section 11111 illustrates a procedure to derive the response to perturbation based on a master equation with a time- 
dependent rate. We also write, using a discrete time picture, the exact solution of this stochastic master equation, as 
a time dependent and stochastic signal. We discuss how to evaluate the statistical averages of this signal in Sections 
IIVI and Ivl In Section llVl we evaluate this statistical average by making the phenomenological assumption that the 
perturbation-induced bias is determined by the tossing of an unfair coin when a collision occurs. In Section Ivl we 
discuss how to convert, in the weak perturbation limit, the perturbed waiting time into a bias. Section IVII illustrates 
the system's response in the time asymptotic limit, in both the phenomenological and dynamic case. In Section IVIII 
we summarize the main results of this paper. The appendixes A and B illustrate the algebra behind the main results 
of Section IVII in the phenomenological and the dynamic case, respectively. 



II. DYNAMIC MODEL 



In this section we discuss how to go beyond the phenomenological approach, by means of a dynamic model that 
has been used in earlier work (see [jj, Il0| and references therein) with different motivations: This dynamic model 
serves the purpose of emphasizing the limits of the Liouville approach to super-diffusion |9j as well as the purpose 
of providing an efficient dynamic algorithm for the non-Poisson renewal condition [T(j- Here we use this model to 
simulate the effect of an external perturbation on a non-Poisson dynamic process, with the following procedure. Let 
us consider the equation of motion: 

V = a(t)y z , (6) 

with z > 1, and a(t) <C 1. As pointed out in Section [I] in this article we adopt a discrete time picture, namely, 

y t +i =y t + a(t)(y t ) z . (7) 

The condition a(t) <C 1 makes it possible for us to convert the discrete into a continuous time picture, insofar as 
the event occurrence involves f > 1, and so an integer virtually indistinguishable from a continuous time. Eq. © 
describes a particle moving within the interval I = (0, 1]: When the particle reaches the border, y = 1, it is injected 
back to new initial conditions within the interval I with the same probability, p(y) = 1. The event of the particle 
reaching the border y = 1 coincides with the crucial collision event that will be discussed in the next section. 

It is important to point out that we have in mind as an example of non-Poisson two-state system the blinking 
quantum dots |20j |. These physical systems are nano-crystals that under the influence of a radiation field produce 
intermittent fluorescence, which is equivalent to a random transition from a "light on" state, with the system emitting 
fluorescent light, to a "light off" state, with no light emission. The waiting time distribution in these two states is 
not exponential, and the "light on" distribution is different from the "light off" distribution. In this paper we are 
interested in establishing the bias created by an external perturbation. Thus, we make the assumption that these two 
distributions, in the absence of perturbation, are identical. In the case when a(t) is time dependent as an effect of an 
external perturbation, its time evolution may depend on whether the system is in the state "on" or "off". Thus, it is 
convenient to adopt the following notation 

a±(t)=a [l + eP±(*)] J (8) 

with the subscripts + and — referring to the states |1 > and \2 >, respectively. The parameter e denotes the 
perturbation strength. The collision event triggers a fair coin tossing, which has the effect of fixing the sign of the 
ensuing state. 

It is worthwhile to point out that the dependence of a(t) on the state sign makes it convenient to distinguish the 
two events. We reserve the term collision to the variable y reaching the border y = 1. We shall explicitly refer to 
the other event as coin tossing event. We shall see that the effect of perturbation (e ^ 0) is to make the collisional 
depend on the coin-tossing event. 
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In the case when a(t) — ceo, namely, when a±(t) is time independent, this model coincides with the generator of 
inverse power law used in Ref. 01 ■ The model used in Ref. is renewal, because the initial condition of y(t) 
after each back injection is selected randomly. Consequently, the ensuing time of sojourn within the interval / does 
not have any memory of the earlier sojourn times. When a is time independent, the model yields the waiting time 
distribution of Eq. (JIJ, with 



T = 



fj, 



and 



z- 1 



(9) 



(10) 



This indicates that perturbing a, while keeping fixed the value of z, is equivalent to assuming that the external 
perturbation changes the parameter T only. 

Let us use Eq. ijBJ) to determine the time r it takes the system, moving from the condition y(t') at time t', to reach 
the border. With an easy algebra, we find 



y(t') 



1 + (*-!) 



T+t' 



a±{t")dt" 



-1/(2-1) 



(11) 



This expression allows us to determine the conditional probability ip±( T \t')i which is the probability of meeting a 
collision at time t' + r given that two events, the collision and the coin-tossing event, occur at time t' . With the 
assumption of uniform back injection, this conditional probability is given by 



4>±(T\t')dT = p[y(t')]dy(t'), 



(12) 



with p[y(t')\ = 1. Note the adoption of a continuous time picture for these calculations is made legitimate by the 
condition 1 <C dr <C T, which is, in turn, a consequence of a(t) <C 1. Thus, we obtain 



^±(r\t') = \dy(t')/dT\ 



a±(r + t') 



The waiting time distribution ip±(r\t f ) satisfies the normalization condition 

/ d,T1p±{T\t') = 1. 

Note that the function ip(r\t') can also be written as follows 

ip±{r\t) = ip±{T\t'), 

with 

t = t' + r. 



(13) 



(14) 



(15) 



(16) 



We note again that the time when the earlier jump occurred, a coin was tossed to decide whether the system had to 
be located in the state |1 >, sign +, or |2 >, sign — . Therefore ip(r\t) is the probability that a jump occurs at time t 
under the condition that the earlier jump occurred at t — r and correspondingly a sign choice was made. 
Using Eq. JSJ) we write Eq. I|13f) in the more convenient form 



a [l + eF±(T + f) 



Let us make some remarks on the physics behind the theoretical prediction of Eq. I|17fl . We note that 

t' = Ti + ... + T„_! 

and 



(17) 



(18) 



T = T n , 



(19) 
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with n — 1 denoting an arbitrary positive integer number. The system is prepared at t = and Eq. (|17(l expresses 
the probability of getting a sojourn time of length r, given the condition that an earlier jump occurred at a time t', 
as the last of a sequel of n — 1 events. 

In the special case of a harmonic perturbation, 



Eq. (jT7|) generates 



F±(t) = ±cos(wt), (20) 



V ' ^ T {l + r/T±^[sin(^(t' + r))-sin(^')]} 



Eq. I|21|) plays a fundamental role for the dynamic approach of this paper. By expanding it up to the first order in 
the perturbation parameter e, we get 

i>±(T\t')tti>(T){l±ecos[u;(t' + T)}}±e^+i{T,t'). (22) 

This first-order approximation maintains the normalization condition, thanks to the second term on the right-hand 
side of Ea. l|22|l . whose explicit expression is of no interest for the discussion of this section. It is enough to say that 
for t — > oo, this second term is proportional to l/r A1+1 . This suggests that the significant contribution to the linear 
response to perturbation depends only on the first term on the right-hand side of Eq. (|22|l . thereby leading to 

i>±(r\t') » ip(r){l ± ecos[Lu{t' + t)}}. (23) 

The approximated expression of Eq. Il2.'{[| suggests a way to convert the exact conditional probability of Eq. I|21|) 
into the joint probability of event occurrence and sign drawing. To properly discuss this issue, let us denote by A t 
and B t the coin tossing and collision events, respectively, occurring at time t. Using the conventional probabilistic 
notations we express i/j±(T\t') of Eq. (|21(l in the following form 

rl>±(T\t/)=P(Bt\At>,Bt>). (24) 

In fact, as earlier pointed out, ip±(r\t') is the probability that a collision occurs at t (B t ), given that fact that the 
earlier collision and the sign choice of the quiescent region occur at t' (A# and Bf). The probability that a collision 
in a state chosen with the coin tossing triggered by the earlier collision (A t >) occur at time t, under the condition that 
the earlier collision occurs at time t', P(A t ' , B t \B t >), fulfills the following equation 

P{A t ,,B t \B t ,)=P{B t \A v ,B t ,)P{A v \B t ,). (25) 

On the other hand, the use of a fair coin yields 



Thus, using Eq. 124(1 . we obtain 



P(A t ,\B t ,)=P{A t ,) = \. (26) 



P(A t ,,B t \B t ,) = ^ ± (T\t'). (27) 



This is an exact result, and the adoption of the dynamic approach should keep us as close as possible to it. 

As mentioned in Section the phenomenological approach to the response of a walker to perturbation rests on 
assuming that the choice between jumping to the left and jumping to the right is done with a biased coin. In the case 
discussed in this article, this corresponds to selecting the state, either |1 > or |2 >, at the beginning of a quiescent 
region. Thus, with the phenomenological approach P(B t \A t >, B t >) = P{B t ) insofar the occurrence of a collision at 
time t does not have any dependence on the earlier sign choice, which does not affect the quiescent region time 
duration. Note that the collision occurring at time t is the first one after the collision at time t'. Thus, due to the 
renewal character of the process, P(B t ) is proportional to ip(t — t'). On the other hand, given the fact that the 
phenomenological approach rests on the sign choice with an unfair coin, P(A t >\Bf) is not independent of Bf. Thus, 
we replace Eq. (12*5)) with 

P(Af,B t \Bf) = P{B t )P{Af\B t ,). (28) 

This yields 

P(A t ,,B t \Bf)=ip(r)p ± (t'), (29) 
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where 

P±(t) = , (30) 

which indicates that we identify P(A t '\B t i) with the unfair-coin tossing event. In other words, the sign choice is done 
at the beginning of the quiescent region, tossing an unfair coin. Note the sign inversion that has the effect favoring 
the state 2 > (|1 >) in the same way as the collision occurring in the state |1 > (|2 >) earlier than in the state 
2 > (|1 >). In other words, as an effect of the approximation of setting e = in the denominator of the term in the 
right-hand side of Eq. Q21[l. with the phenomenological approach we assume that the conditional probability ip(T\t') 
of Eq. (|21|l yields the joint probability 4 ! ( T )p±(t / ), which is the probability that at time t' we assign the sign + or 
— to a quiescent region of time duration r. The two events, the drawing of time r and the sign selection, are totally 
independent. As a consequence, there is no correlation between sign and time duration, a fact in a striking conflict 
with the well known fact that stochastic resonance can establish a strong correlation between sign and time duration. 
The reader can consult Ref. |21| for a case of strong departure of the perturbed waiting time distribution from the 
unperturbed form. 

Let us consider now the case where the biased sign choice is made at the end of the quiescent region, namely, the 
case 



P(A t ,B t \B t ,)=i,(T) P± (t), (31) 

with 

1 ± ecos(wt) .„„. 

P±(t) = g"-^. (32) 

Note that in this case there is no sign inversion. This choice is equivalent to 

P(A t ,B t \B t ,)=xP ± (T\t')p±(t), (33) 



where ip±i T \t') is given by Eq. (|23|l and 



P±{t) = \. (34) 



Thus the adoption of this choice allows us to maintain the promise of remaining as close as possible to the exact dynamic 
prescription of Eq. J57J). This can be realized in two distinct but equivalent ways. With the first prescription, we toss 
a fair coin either at beginning or at the end of the quiescent region, whose time duration, as a perturbation effect, 
depends on whether the system is in the state |1 > or |2 >. With the second prescription we select the sign of the 
quiescent region, whose time duration is unperturbed, tossing a unfair coin at the end rather than at the beginning 
of the quiescent region. With this prescription we depart from the phenomenological approach, where the sign choice 
is done, tossing an unfair coin, at beginning of the quiescent region. 

The dynamic and the phenomenological approach, which coincide in the Poisson case, lead to different predictions 
in the non-Poisson case, and the departure of one prediction from the other becomes significant especially when the 
condition fi < 2 applies. This is a consequence of the fact that the adoption of a biased coin, tossed at the beginning 
rather than at the end of the quiescent region, does not reflect quite properly the dynamic origin of the response to 
perturbation. If the phenomenological model is thought of as an approximation to the dynamic model, it is evident 
that it has to be limited not only to the case of extremely weak perturbations, as the dynamic approach, but also to 
cases where the perturbation-induced waiting time reordering can be totally neglected, so as to fulfill Eq. JSJ). 

III. STOCHASTIC LIOUVILLE EQUATION 

Let us consider the ordinary master equation 

p(t + l)-p(t) = ~Kp(t), (35) 



where 



K=( \ / ) C'Xi! 
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and p(t) is a vector, whose components pi(t) and p-2.it), with p\(t) +pi{t) — 1, denote the probability for the system 
to be found in the corresponding states |1 > and |2 >, respectively. We are considering the case of a discrete time 
t. Thus, Eq. l|3*5f is nothing but the probabilistic picture corresponding to a fair coin-tossing prescription. We can 
imagine a particle located, for instance, in the state 1 >. If the fair coin tossing yields head, the particle remains in 
this state. If the fair coin tossing yields tail, the particle jumps to the state 2 >. If the next coin tossing yields tail 
again, the particle remains in the state |2 >. Otherwise it jumps back to the state |1 >. We define the coin-tossing 
event as collision. The time distance between a collision and the next is equal to 1, in this case. Thus, collisions are 
not a source of disorder, which, in this case depends only on the coin tossing. 

It is more realistic to assume that the time distance between one collision and the next is not fixed and fluctuates. 
The leading idea of Continuous Time Random Walk (CTRW) can be adapted to this simple process by assuming 
that the time distance between one collision and the next is not constant, and it is characterized by a probability 
distribution ?/>(r). In this paper we develop a theory for the non trivial case where the key function ip{ T ) nas a 
non-exponential form. We shall make explicit calculations for the special case where - 0( r ) has the form of Eq. 
Note that the theory developed in this paper makes it possible for us to address also the case [i < 2. In this case, the 
mean waiting time is infinite, and the system is found to be in a state of perennial aging, insofar as the stationary 
condition is not admitted in this case |15|, [la, 113, LL8J • 

The theory of this paper rests on the connection between the CTRW and the Stochastic Liouville Equation (SLE) 
of Kubo [23L I24I |25| . Let us devote some attention to this crucial aspect. As a consequence of the time fluctuations 
between one collision and the next, we write Eq. ll-iol) under the following form 

^~=*W (37) 

where f(t) is a two-dimensional vector with components fi(t), f2(t) and fi(t) + fift) = 1. The fluctuating rate ro(t) 
always vanishes but in the correspondence of a collision, where it gets the value of 1. The functions f\(t) and /2(f) 
denote the probability for the particle to be in the corresponding states |1 > and |2 >, as in the earlier case of constant 
time distance between two consecutive collisions. The reason for this notation change is that, as in the SLE by Kubo 
|25|, fi(t) and f2(t) are fluctuating probabilities, rather than ordinary probabilities: In the original work of Kubo 
|25j the fluctuating probability f(x,t) is turned into an ordinary probability, obeying the Fokker-Planck equation, by 
averaging on the stochastic process. 

It is important to stress that our approach is based on a discrete time picture, with the integration time step dt = 1 . 
Therefore, we define 

^ = f(i + l)-f(t). (38) 
Consequently, when a rare collision occurs, Eq. i(3T|) produces the same effect as Eq. (|3*5|) . Let us define 

m = Mt) - h(t). (39) 

Thus, Eq. (J23 becomes 

E(t) = E(t-l)-r (t-l)E(t-l). (40) 

Note that 

S(t+1) = (41) 

and 

E(f + l) = E(f), (42) 

according to whether at t a collision does or does not occur, respectively. Thus, a collision occurring at t turns E(t) 
into the vanishing value, while, with no collision, the value E(t) is transmitted to the next time step. 

The earlier approach is a convenient way to describe the relaxation process in the unperturbed case. To address 
the problem of the response to external perturbation we are forced to adopt the following generalization: 



f(f + l)_f(f) = _± 



r+(t) -r_(t) 
-r+(t) r_(t) 



f(t). (43) 
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With r±(t) we denote a stochastic process that is almost always vanishing, except at the moment of a collision. If the 
collision occurs at the same time for both the |1 > and \2 > state, phenomenological model, the bias is determined 
by the assumption of Eq. (|32|) , which can be interpreted as generating the state dependent rate 



r±(t) = r (t)[l±ecos(wt)], 



(44) 



where ro(t) is a stochastic function of time, which is almost always a vanishing quantity, except in the rare case of a 
collision, when it gets the value of 1. The time distance between two non-vanishing values of ro(t) is randomly drawn 
from the distribution iP(t) of Eq. (JTJ. According to the dynamic models, the collisions in the states |1 > and \2 > 
occur at different times, and, consequently, r+(i) and r_ (t) are equal to 1 when the collision occurs. Let us imagine, 
for example, that the collision occurs in the state |1 >. We see that Eq. 1431) corresponds to shifting half particles 
from the state |1 > to the state \2 >. The opposite effect is generated by a collision occurring in the state |1 >. 
In conclusion, to take into account the effect of an external perturbation we have to turn Eq. (|40|l into 



E(t) = 



1 - 



S(t-1) 



E(t-l)- 



D(t-1) 



where 



and 



D(t)=r+(t)-r-(t) 



S(f) = r+(t)+r_(f). 



The solution of Eq. (jHSJ is 



where 



if t = t' + 1, and 



£(i) 



t'=Q 



Q(t,t') = 



1 - 



Q(t,t') = l, 



S(t' + 1) 



S(t-1) 



if t > t' + 1. Using Eq. we write Eq. igSJ as follows 



t-i 



m = - £ 



Q(t,t') 



t'=o 



t-i 



t'=o 



(*') 



Q(t,t'). 



(45) 



(46) 



(47) 



(48) 



(49) 



(50) 



(51) 



It is important to point out that we prepare the system. This means that r + (0) = r_(0) = 1. Thus E(l) = 0. 
In the continuous time limit this is indistinguishable from E(0) = 0, which is equivalent to assuming that there is 
no bias on the initial condition. We note that for t > 1, the value E(t) is determined by the events occurring at 
times t' meeting the condition t — 1 > t' > 1. It is important to notice that Q(t,t') depends only the occurrence or 
non-occurrence of events at times larger than t' and smaller than t. 

The key issue to discuss in the next two sections is as to how to evaluate 



II(t) =< £(i) > 



Using Eq. (|51|l . we write 



n(t) 



i t_1 

< E(t) >= -- [< r+(t')0(«,0 > - < r-(t')Q(t,t') >} 



(52) 



(53) 



t'=o 



The quantities < r±(t')Q(t,t') > are evaluated over infinitely many trajectories. All these trajectories begin with 
the beginning of a quiescent region, half of them being initially located in the state |1 > and half in the state \2 >. 
Thus, the quantities < r±(t')Q(t,t') > can be interpreted as being the probabilities that the sequences of events 
corresponding to the definitions of r±(t')Q(t,t') really occur. We shall see that these probabilities can be expressed 
in terms of the quantities P(i), ^(t), ip(t, t'), all of them being derived from the distribution ip(t) of Eq. (JTJ. 
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IV. PHENOMENOLOGICAL APPROACH 



In this case, we use the prescription of Eq. ||4"4T . Thus Eq. yields 

t-l 

11(f) = -e ^ < r o{t')Q(t, t') > cos(wt'). (54) 

t'=0 

The time distance between two consecutive events, namely, two occurrence times where ro(t) is equal to 1, has the 
distribution ip(r), because we are assuming that the external perturbation does not affect the time distribution of 
events but only the rates r+(t) and r~(t) when an event occurs. 

Let us now focus our attention on the function < Q(t,t') >. We want to prove 

< Q(t,t') >=V(t,t') = V(t-t'). (55) 

Note that we rest on the probability that at t' an event occurs. If it did not, ro(t') = 0, and there would be no 
contribution to the sum of Eq. (|54|l . The function Q(t,t'), according to the definition of Eqs. I|49l) and i|50|) . is the 
product of non vanishing factors with the value of 1, if no event occurs in between t' and t. In fact, if no event occurs 
at t' + 1, S(t' + 1) = and 1 — S(t' + 1) = 1. The same argument applies to the following times until time t — 1. At this 
time no event occurs either, thereby making S(t — 1) = and 1 — S(t — l)/2 = 1. Thus, if no event occurs in between 
t' and t, as earlier stated, Q(t, t') — 1. Let us assume now that an event occurs at t. As a consequence S(t) = 2. Thus, 
at the next step we get 1 — S(t)/2 — 0, thereby making Q(t + l,t') = 0. Thus Q(t,t') is equal to 1 and it vanishes 
when we make the change Q(t,t') — > Q(t + l,t'). On the other hand, at t 1 an event occurs. Thus, < Q(t,t') >, the 
probability that Q(t,t') does not vanish, is equivalent to the probability that we do not find any event moving from 
t' , when an event occurs, up to time t. Consequently, not only < Q(t,t') > is a survival probability: it is a brand 
new survival probability, namely evaluated by setting the beginning of the observation at a distance At — 1 from 
the occurrence of the last event. The measure of this distance becomes virtually zero in the continuous time limit, 
thereby explaining why we obtain Eq. I|55|) . 
In the continuous time limit we write 

U(t) = -e [ dt'x(t,t')cos(ut'), (56) 



which is the typical linear response function structure 0,0 • In literature, the function t') is called susceptibility. 
Its explicit expression in the phenomcnological case here under discussion is: 

X(t,0 = *(*-*W), (57) 

where 



P(*) = £>n(*), (58) 



with ip n {t) denoting the probability of occurrence, at time t, of the last of a sequence of n collisions. Note that the 
sum is done from n — 1 to oo, as a consequence of the fact that n = corresponds to the preparation event located 
at t — 0. This event does not contribute to the bias, insofar as we assume that at t — there are as many particles 
in the state |1 > as in the state |2 >. Thus, P(t) is the probability of event occurrence at t, given the condition that 
the system is prepared at t = and the event occurrence is observed at time t > 0. The Laplace transform of P(t), 
P(u), is given by 

P(u) = Mu) = - 1 = t^ttt- (59) 

n=l n=0 1 - W) 

where if>(u) is the Laplace transform of - 0( T ) — V'i( T )- I n the case \x < 2, the Laplace transform of iP(t) of Eq. Q for 
u — > has the form 

tp{u) = 1 -T(im- l){uT) t "- 1 . (60) 
By plugging this form into Eq. (|59|l and inverse Laplace transforming the result, we get that for t — » oo 

P(t) cx (61) 
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Thus, the rate of event occurrence decreases as a function of time, an evident aging property that is taken into account 
by the non-stationary susceptibility of Eq. J57J). 

The meaning of the susceptibility t') stemming from the phenomenological condition, is evident. It is propor- 
tional to the probability of event occurrence at t' by the probability that no further event occurs up to time t. Thus 
the bias H(t) collects all the collisions occurring at earlier times and establishing the current P(t') cos(wt'). Note that 



tl( u ) = -e^^Re | $2±M \ . (62) 
u [[I- l/j( U + iuj)} J 

In Appendix A we illustrate the procedure necessary to inverse Laplace transform this formula in the asymptotic 
condition u — » 0, so as to evaluate H(t) for t — ► oo in this case. 



V. DYNAMIC APPROACH 



As noticed in Section [HI the dynamic approach derives the system's bias from the fact that, as an effect of 
perturbation, the events in the state |1 > and the events in the state \2 > occur at different times. Let us imagine 
that two independent trajectories enter the states |1 > and \2 >, respectively, at the same time t = 0. Thus, at t = 
two simultaneous collisions occur in both states. Let us denote by w+(t') and W-(t') the probabilities that a further 
collision occurs at time t' in the states |1 > and 2 >, respectively. Using Eqs. (|33|l . (|34|l and l|23|) . we predict that 

w+(t') = l + ecos(uit') 

w-(f) 1 - ecos(wt')' 1 ' 

As a consequence of this equality, a collision occurs in the state |1 > earlier (later) than in the state |2 >, if cos(wi') > 
(cos(cjt') < 0). This has the effect of decreasing (increasing) the number of particles in |1 > and increasing (decreasing) 
the number of particles in |2 > by the same amount. 

To convert the occurrence of collisions in the two states at different times into a linear response to the external 
perturbation, let us proceed as follows. Let us set r_(t) = r + (t) = 1 and let us write Eq. (|5*TJl under the following 
form 

=(*) = - E r -^^Q(t, o + E r -^^Q(t, o. (64) 

t'=0 t'=0 

The first term on the right-hand side of this equation can be interpreted as the occurrence of two events, a first event 
occurring in the state 1 > at time t' and a second event occurring in the state 2 > at a later time t > t'. In fact, 
if r + (t') = there would be no contribution to this term. This term would also vanish if r_(t) = 0. But this would 
violate the condition that this term coincides with the first term of Eq. I|51fl . With the same argument we prove that 
the second term on the right-hand side of this equation can be interpreted as corresponding to a first event occurring 
in the state |2 > at time t', and to a second event occurring in the state |1 > at a later time t > t' . We have now 
to move to the crucial step of making a statistical average. It is evident that in the first term on the right-hand side 
of Eq. (|64l) the existence of the terms of the sum over t' depends on the probability that r_(t) = 1 and r + (t') = 1, 
and that, analogously, in the second term the existence of the terms of the sum over t' depends on the probability 
that r + (t) = 1 and r_(t') = 1. Let us make the assumption that no event occurs in between t' and t, thereby setting 
Q(t, t') = 1. The average of E(t) of Eq. @U yields 



n(t) = -\ 



dt' 'M (+ _ + _ ) (t') - / dt'M(—> +) (t') 



(65) 



where Mt±^,^(t') denotes the probability that at time t' an event occurs in the state |1 > (|2 >) and another event 
occurs at time t in the state |2 > (|1 >). Note that in the unperturbed case 

M (+ ^_ ) =M ( _^ +) =^(i,0- (66) 

This is so because the occurrence of an event in the state |1 > (|2 >), at time t', signals the beginning of observation 
in the state |2 > (|1 >). In both states an event occur at time t and observation begins at time t' . On the other hand, 
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to properly take into account the effect of perturbation, we use Eq. I|63|l . Thus, due to the correlation between event 
occurrence and quiescent region sign, we obtain 

M i+ ^-)(t') = w+(f) = [1 + ecosMQ] 
M(__ +) (#) w-(t') [1 -ecos(wt')]' 

To make Eq. I|67[l compatible with Eq. (|66l) when e = 0, we set 

M (+ ^_) =< r + (t)r-(t')Q{t,t') >=< r {t)f(t,t') > [1 + e cos(cjt')] = ip{t, f)[l + ecos(urf')] (68) 

and 

M ( _^ +) =< r-(t)r + {t')Q(t,t') >=< r (t)^{t,t r ) > [1 - ecos(cjt')] = 4>(t,l/)[l - ecos(wi')]- (69) 

It is important to observe that these two equations correspond to the prescription of Eqs. (|33|l and (|34|l . involving 
Eq. I123H as well. As pointed out in Section |3 this prescription can be realized in two different ways, one resting 
on the unfair coin tossing occurring at the end of the quiescent region, and the other based on the use of a fair coin 
tossing and of a perturbed waiting distribution. Here we are using the second way of realizing this prescription. 
Using Eqs.JHHJ an d l|69|l . we rewrite Eq. (|65|l as follows 

H(t) = -/ dt' ^ < " t ' t '\ l + ecos(ujt')}+ f dt' ^^ [1 - ecog(hrt')] = -e / dt'x{t, t') cos(a;t'), (70) 
Jo 2 Jo 2 Jo 

which yields the same form as Eq. 1)56(1 with the susceptibility of Eq. (|70|l replaced by 

X (t,t') = iP(t,t'). (71) 

This result is based on the assumption that no event occurs in between t' and t. On the other hand, this approximation 
can be expressed in a different way. We made the assumption that the bias at time t depends only the delay between 
the event on |1 > (|2 >), occurring at time t and the event occurring in the other state |2 > (|1 >) at time t' . This is 
equivalent to assuming that at time t' there is no bias. If we relax this assumption and consider also the possibility 
that at time t' there is a bias, this is certainly proportional at least to e. Thus, this correction to the approximation 
adopted to get Eq. l(7U|) would be proportional to e 2 . In conclusion, the assumption yielding Eq. (j7T|l is a linear 
response assumption. 

It is interesting to notice that in the Poisson case the prescription of Eq. H7Ufl coincides with the prescription of Eq. 
(156)1 . Notice that the theoretical result of Eq. 1)70)1 coincides with the prescription of Ref. [l^- To make this paper 
as self-contained as possible, let us review the calculation done in 0] to deal with the case \i < 2, which makes aging 
a perennial condition of renewal systems. It is known 01 ( see a ls° R-ef- G3) that the exact expression for the aged 
ip(t,t')is 

tP(t,t') = m+ I P(rMt-T)dT, (72) 



o 



where P(t) is the time-dependent rate of event occurrence of Eqs. I|58|l and (|61fl . the latter of which makes evident 
the perennial aging condition of the case [i < 2 here under discussion, and evidently properly considered by Eq. 1)710. 

It is straightforward to evaluate fl(n), the Laplace transform of of II(t) of Eq. JTHJ). This is so thanks to the time 
convolved nature of this equation and to Eqs. I|59|) and (|61|) . We obtain the following expression 

fl(u) = -eRe(E(u)), (73) 

where E(u) is the Laplace transform of 



E(t)= / dt'tp(t,t')exp(-iojt'). (74) 
Jo 

After some algebra, detailed in Appendix B, we find 



lJj(lL + ILO) — 1p{u) 



E(u) = - ; T 2 - ■ (75) 



1 — ijj(u + iuj) 



In Appendix B we give also details on the procedure adopted to inverse Laplace transform this important equation. 
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VI. ASYMPTOTIC BEHAVIOR 



We are now in a position to derive the central result of this paper. By means of the calculations illustrated in 
Appendix A we show that the adoption of the phenomenological theory yields the asymptotic condition 

n t =ni t +n 2 t we ; -— - — T sm— + — '-. (76) 

w w w \T + t) r(2 - ^{uty- 1 2 r(/x-i) {uot) 2 -^ y ' 

In the dynamical case, as shown in Appendix B, we recover the result of Ref. |l8j: 

e cos (^ + tot) , x 

n (t) « — '-. 77 

Note that in both cases we do not adopt the restriction to the condition uiT <C 1, a constraint used in Ref. [lg. We 
have also to point out that we are considering the time asymptotic condition cot ^> 1. 

For 1 < fj, < | the first two terms in the right-hand side of Eq. (|76[1 are the dominant ones, while for | < \i < 2 
the dominant is the third term. Thus, when /i is very close to 1, the phenomenological theory yields a result with 
no apparent resonance effects, namely, a significant departure from the prediction of the dynamic approach. On the 
other hand, we note that tuning the external perturbation to the frequency 



UJ = U! c — 



sin Hi " ~ 



r(2- M ) 



^ (78) 



has the interesting effect, in the limit t — > oo, of making the first and second term on the right-hand side of this 
equation cancel with each other, thereby making Eq. I|76[) become identical to Eq. (|77|l . and the phenomenological 
identical to the dynamic prediction. At the moment of writing this paper, the physical meaning of this critical 
frequency is not quite clear. 

It is of some interest to discuss the step-like perturbation 

F±(t) = 0(t). (79) 

In both phenomenological and dynamic case the response to a step-like perturbation yields a finite asymptotic value 
for n(£). In the dynamic case this asymptotic value, though, may become much smaller than in the phenomenological 
case. In fact, the phenomenological case yields, see Appendix A, 

n(oo) = -e (80) 

and the dynamical case produces, see Appendix B, 

E(oo) = -e(M-l), (81) 

which vanishes for fj, — » 1. 



VII. CONCLUDING REMARKS 

This paper affords a detailed account of the dynamical theory proposed in Ref. |l8j |. Furthermore, the time 
asymptotic behavior of the response is evaluated with no restriction to the condition ujT <C 1, adopted in Ref. 18]. 
It is interesting to notice that with the dynamic approach a form of resonant response is admitted, even if it slowly 
dies out. With the phenomenological model the non-resonant response may dominate, preventing us from observing 
any form of resonance. The response to a step-like perturbation, in the time asymptotic limit yields a finite sig nal. 
Note that the time asymptotic value of the response to a step-like perturbation, in the case studied in Ref. |l5l Utij 
is a vanishing current, in an apparent conflict with the prediction of this paper, where both phenomenological and 
dynamic approach yield a non- vanishing current, although of different intensity. As pointed out in Section [I] our 
approach refers to the super-diffusion case, thereby explaining the departure from the results of Ref. |!5lll6|. 

In conclusion, we think that this paper, on the one hand, affords the possibility of establishing whether a given 
system obeys the dynamic or the phenomenological model. Should further investigation, for instance in the field of 
blinking quantum dots 20], prove that the dynamic model is more plausible than the phenomenological model, the 
procedure of this paper may be adopted to derive the phenomenological model from dynamics, and to assess under 
which conditions the phenomenological model can be trusted. It seems that the approximation adopted to derive 
the phenomenological condition from dynamics becomes negligible with fi — > 2. On the other hand, we think that a 
challenging problem is to overcome the weak perturbation approach necessary to make theoretical predictions with 
the dynamic model. 
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APPENDIX A 

This appendix is devoted to inverse Laplace transforming Eq. i|62[) that we rewrite here for the reader's convenience: 



II (u) = -eRe 



1 — t/j(u) 1p(u + IU)) 



u 1 - %j)( U + uj) 

We express the Laplace transform of ip(u) by means of the exact expression |26|,|2j 



= r(l - H) {Tuf- 1 (// - 1) (e uT Ef_,) 



where 



pi _ \ x 

7 ~ ^ r(n + 1 



n=0 v ,; 



Using Eq. ljA-21) . let us write Eq. I|A-1|) as: 



n (u) = -eRe 



1 - ?/>(«) $(tt + vJ) - 1 + 1 



1 - #0 R 

e cfte 



1 1 - 
U 1 - + w*>) 



(A-l) 



(A-2) 



(A-3) 



(A-4) 



u 1 - ^>(it + «w) 

Let us set f[ (u) — tl± (u) + fto (u) so that LTi (u) — e 1- ^^ and IT (u) = —eRe — • — ; . It is straightforward 
to inverse Laplace transform ITi(w), which yields IL (i) — e ( jq^-J ■ The inverse Laplace transform of IT2 (u) is 

more laborious. First of all, we introduce the notation IT2 (u) — Re lie i u ) 
right side of Eq. (|A-4|) as 



and we rewrite the second term in the 



tic (u) = - e -{u + WW 

u {u + iuy-f* -r(i -ii)T»- i {n- 1) e («+»)T_£ "^) T 



(A-5) 



{(it + iw) 1 "' 1 - r(l - ^T^ 1 ^ - 1) \e^ u+l ^ T - ^«+~) T ] J n c (u) = _£( u + kj)!-^ [1 - 

To proceed with our calculations, we have to prove first the important formula 

00 

t a 



t + T 



-e ~ ut dt = 



sin 7ra 



(A-6) 



(A-7) 



Let us replace the exponential exp(— ut) on the right-hand side of Eq. (|A-7|) with its Laplace transform, with respect 
to u, 



(A- 



P + t 



We are therefore led to the integral |28g 

,00 ta 1 



-dt 



q t + T t + p — T p — T J sin(7ra) 



(A-9) 



According to the Efros theorem |29|, the inverse Laplace transform of Eq. I|A-9(1 . from p to u, is the Laplace transform 
of t a /(t + T). Using Ref. [27j we derive from this procedure the right-hand side of Eq. i|A-7|) . This proves the equality 
ofEq. 
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Going back from Eq. (|A-7|) to the time dominion we have: 



1 



T(p-l) 



t'^e-^'Uc (t-t')dt' 



1 



rou-i) J f + T 



e~ lult 'li c (t-t')dt' = git), 



(A-10) 



where g{t) is the inverse Laplace transform of the right side of Eq. l|A-6j) . After a little algebra: 



T 



T(p-l)J t' + T 



t HL- 



. e Mt-*)Uc (t - t') dt' = e luJt g{t). 



(A-ll) 



The Laplace transform of the left side of Eq. IjA-llf) is k(u)Ilc(u — kj), where k(u) is the Laplace transform of the 
kernel. As to u — > 0, we write k(u) = fco + w(u), where 



T 



r(M-i) J f + T 





t lfL- 



dt' 



7T T^- 1 
sin7r(/i) T(fi — 1) 



(A-12) 



We do not write the explicit expression of w(u) since w(u) — > for u — > 0. The Laplace transform of the right-hand 
side of Eq. (| A- 1 1|> is g(u — iuj); dropping off the ilo translation, we study directly the equation for TLc(u) and g(u), 
so that we end up with: 



k flc(u) = (u + iu;) 1 -^ l-4>(u) 

u 



(A-13) 



Taking the limit u — > we arrive at: 



n c («) = r(i - ix)t^ 1 ^- — — — ( -^—Y 1 = e m - ^r*- 1 ^— ^"" 2 

ko u \u + %lu J ' ko \ u + ilo 



1 



(A-14) 



It is possible to express Eq. (|A-14ll in terms of a fractional derivative. We use the definition of fractional derivative 
of Ref. [13, which yields 



c- 1 



u-2 



1 



U + ILO 



fi-2 



^— 'lLdtj.fi — 2 

r(M - 1) 



and, consequently, 

c- 1 



n— 



r(M-i) 



Using the expression for fco of Eq. I|A-12I) . after a few calculations we obtain: 



n=0 



n /r(n+l) 



-eF(2 - /i, l,iwi)e~ 



which can be expressed for t — > oo as [26 

n c (t) 



r(2- M ) 

Finally, using the real part of ilc(i), we have 



r( M -i 



(A-15) 



(A-16) 



(A-17) 



(A-18) 



n (t) = ni (t) + n 2 (t) » e 



T 



T + t 



■ sm 



7T/I 



e cos 



F(2-/z)(o;t)' i - 1 2 r(/i-l) (wt) 2 "' 1 



(A-19) 
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which is the key result of Eq. (|76|l . 

The response to the step-like perturbation of Eq. (|79|l is obtained by the earlier results, by setting u> = 0. In this 



n(u) = -e 



tp{u) 



and 



H(t) = -e 



1 - 



T 
T + i 



(A-20) 



(A-21) 



which proves Eq. (|50|) . 

We stress that in this Appendix we did not set no constraint on T and ui, as done in the earlier paper of Ref. [l8| . 
where the condition wT <C 1 was adopted. Neither the calculations of Appendix B, made along the same lines of this 
appendix, are forced to obey the condition luT 1 of Ref. 0] . 

APPENDIX B 

In the first subsection of this Appendix we derive Eq. I|75|l . In the second subsection we show the algebraic details to 
inverse Laplace transform this equation in the asymptotic limit, and we find the response to the step- like perturbation 
as well. 



A. Analytical Laplace Transform 

Let us make the Laplace transform of Eq. (|74|l . We apply the method of integration by parts, yielding 



E{u) = I i 
lo 



l dt \ ip{t, t')e- tut dt 1 = 
Jo Jo 



1 1 d( e - iuJt ') 



This means: 



E{u) = - 



i 



'dt 



-iio dt 



dt' 



■j—^ipit^t^dt' . 



Using the definition of ip(t,t'), given by Eq. I|72|l . we get: 

^(t,t)e- iut - i>(t) - / e- iut ' J2Mt'Mt-t')dt' 
Jo n=l 

Using again the definition of Eq. (|72[1 . we have 



E(u) = - / e~ ut dt 
u Jo 



ui 



E(u) = - / e~ ut dt { e 



oo „t 

/ MrH(t~r)dT 

n=l J 



(B-l) 



(B-2) 



(B-3) 



m- e- iut 'J2Mt')i>{t-t')dt'\. (B-4) 



Taking into account that we are dealing with uncorrelated events (ipn{u + i<^>) = (V-K" + iw)) n ) we write the Laplace 
transform of Eq. (|B-4(1 as follows: 



% i * x — > r * i n * ~ x — r * i n A 

E(u) = — < + iuj) + ^P( u + w) i>(u + — tp(u) — + iu) ip(u) 



I 



oc n 

ijj(u + iuj) — V>(it)J + \^>{ u + ' + * w ) — 



(B-5) 
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By summation of the geometric series, we obtain the important result of Eq. J75J, written here again for the sake of 
reader's convenience: 



- i 4>(u + iu) - V>(u) 

E[u) = — • . 

w 1 — ip{u + iui) 



(B-6) 



As far as the step function case of Eq. (|79|l is concerned, we can find the corresponding Laplace transform with 
two different methods. The first rests on evaluating the limiting case of u> — > of Eq. (|B-6|I . Thus, we obtain: 



E(u) 



1 ip(u + iui) — ip{u) 



1 



lim 

ui-rQiLj l-ip( u + iuj) l-ip(u)du 



(B-7) 



As a second method, let us evaluate the Laplace transform of J Q ip(t,t')dt', namely, the response generated by the 
step function, to prove that it yields Eq. (|B-7|I . We write 



e~ ut dt / ip{t,t')dt' = / e~ ut dt\ [t'ip{t,t% - / t 
o Jo Jo I Jo 



E{u) 

By applying a method similar to the one used to derive Eq. (|B-3J| . we obtain 

E(u) 



1 ,dmx)) dt , 



dt' 



(B- 



p oo 




/ e~ ut dt 


tip(t)+t f 


10 


Jo 



ti>{t)+t y^i> n (T)i>{t-T)dT- I ^t'^ n {t')4,{t^t')dt l 
JO n=l JO n=l 



(B-9) 



By making the Laplace transform of all terms on the right-hand side of Eq. (|B-9J) . we finally arrive at: 

run- +^)iE^r = -7T^ d 

n—1 J n—1 



E{u) = ~i>(u) j- 
du du 



^(u)j^(u)' 



tp{u) du 



(B-10) 



which is the expected result. On the basis of this result, we conclude that the response in time to the step function 
perturbation can be derived from the general result of the harmonic case for ui — > 0. Thus, at the end of the next 
subsection we shall use this argument to find the response in time to the step function perturbation. 



B. From the Laplace domain to the ordinary time domain 

At this stage we must make Eq. (|B-6|) suitable for the inverse Laplace transform procedure. Therefore we have to 
find the analytical form of ip( u ) f° r u —>■ 0. Let us use again the same departure point as in Appendix A, namely Eq. 
(|A-2() . which we rewrite here for the sake of reader's convenience 



$(u) = r(l - M ) (Tuf- 1 (fj, 1) (e« T - Ef^) . 
With the help of this important property, we rewrite Eq. (|B-6() under the following form: 

1p(u + lUj) — 1p(u) 



E(u) 



-{u + iuj) 1 -^ 

w (u + iujy-t* -r(i -/x)tm-i( m - !) e (u +l ^)T „ E i^)T 



From this equality we get: 

[{u + ilo) 1 -^ - F(l - - 1) 



For the sake of clarity we rewrite Eq. (|A-7J| 

oo 

/ 



t° 



t + T 



sin na 



-T a (E l ' T -e 1 ' 



(B-ll) 



(B-12) 



e («+«*)T _ s (u+MTj } E(U) = 1 {U + ujf-i* \j^ u + - ^(«)] . (B-13) 



(B-14) 
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Following the same procedure of Eq. (jA-lOfl . going back to the time dominion, we have: 

t t 



r(M-i) 



t' 1 '- 



-2 „ — iujt 



E(t- t') dt' 



1 



r(/i- i) J t' + T 

o 



e~ iut E{t-t')dt' = g(t), 



(B-15) 



where g(t) is the inverse Laplace transform of the right side of Eq. I|B-13|I . After a little algebra, we arrive at: 



T 



V{ji-l)J t' + T 



E(t-t')dt' = e lu,t g(t). 



(B-16) 



The Laplace transform of the left-hand side of Eq. i|B-16|) yields k{u)E{u — iuS), where fc(it) is the Laplace transform 
of t^ 2 / (t + T). In the limiting case u^Owe write k(u) = kg + w(u), with 



fcn = 



T 



■dt' = -- 



T(fj, -1) J t' + T sin7r(/x) - 1) 



(B-17) 



We do not write the explicit expression of w(u), which tends to vanish for u — > 0. The Laplace transform of the 
right-hand side of Eq. I|B-16(1 is g(u — iu>). We do not make any additional approximation, and we study directly the 
equation for E(u) and g(u). After a few calculations we get: 



k E (it) = — (it + id) p i/)(u + uS) — ip(u) 



(B-18) 



Using Eq. (|B-11|I we finally obtain: 



E(u) =r(l-n)T> 



ko w 



Ju+iuj)T _ g(u+iu)T _ 



At-1 



u + iu> 



(B-19) 



We stress that so far only the function A; (it) is ap prox imated while the right side of Eq. I|B-18() is exact and we do 
not set any constraint on T and u, as done in Ref. [lq. From Eq. I|B-19|) . taking the limit for u — ► and w->0we 
obtain: 



II(it) = eRe 



u + iuj 



(B-20) 



This equation was found in Ref. [18j, under the assumption ujT -C 1. Let us invert Eq. I|B-19J) term by term. Using 
Eq. HB-14JI . we obtain: 



sin7r(/x - 1) t^ 1 
ttT^- 1 t + T e 



(B-21) 



and 



It + IUJ 



(e uT - Ef_ t 



r(n-i)T(i-n)j (t-t' + TY 

o 



e- luJt 'dt'. 



(B-22) 



Let us develop the exponential factor of the right-hand term of Eq. 1B-22(1 into a time series, and let us evaluate the 
corresponding integrals. This yields: 



ZT 1 



u + uo 



K T -i) 



v -FU,n + /i-l,n + M,-— , (B-23) 



7r(t + T> ^ (n + u - l)n! 

v y n— 



18 



where F(a, j3, 7, z) is the hypergeometric function. Using the results of Eqs. (|B-21() and (|B-23I . after a little algebra, 
we get for E(t) the following expression: 



E(t) = - 



(/1 — l)t p 1 sin7r/z 



E 



n(t + TY ^ (n + fi — 1) 



(-lU)t) n ( t 

:F I fi,n + fi — l,n + Mi 



t + T 



sinTT/^ 1 - e" 4 "* / i 
7T «wT V T 



n-x 



T 



t + T 
(B-24) 



In the time asymptotic limit t — ► 00, it is possible to evaluate the sum appearing in Eq. l|B-24fl . This yields: 



E(t) = 



t + T 



e" lut {n - 1)F(2 - n,2,iut) + 



sin7r fi 1 — e 41 

7T IwT 



T 



/i-2' 



sin7r/x 1 — e 11 

7T «wT 



t\^ 2 t 



T 



t + T 
(B-25) 

where F(a,j,z) is the confluent hypergeometric. The first term between the square brackets, multiplied by — e, 
coincides with Eq. (34) of Ref. J^. We know |l| that for t -> 00 F(2 - fjt, 2, tut) ~ ^ yielding Eq. (36) of Ref. 
^8|. We now prove that this earlier result can be recovered without setting the condition uiT <C 1 adopted in Ref. 
|18|. In the time asymptotic limit, Eq. (|B-25|I yields: 



E(t) w (/x - 1)^(2 - M> 2 , 



sin7r/i 1 — e u 

7T ZWT 



T 



/i— 3 



(B-26) 



For cj — > the second term of the right side of Eq. (|B-26|) can give a non-negiglible contribution in a certain region of 
time, more precisely when t — ► 00 but still with tot <C 1. In fact, in this case the exponential function in the second 
term of the right-hand side of Eq. (|B-26(I produces a factor of the order t/T. Moving to n(t) we obtain: 



, . „ r . „ „, , t1 sin nu sin u>t ( t 

7T LJ1 \ 1 



fi—3 



(B-27) 



with a dependence on T. For u;i ^> 1, taking into account that /i < 2 and the asymptotic properties of the confluent 
hypergeometric function j2?| as well, we get: 



n(t) « e 



cos + wt) 

r(^-i)(wt) 2 -/*' 



(B-28) 



which coincides with Eq. (36) of Ref. without involving the condition uoT <C 1. Note that Eq. (|B-28|I is the 
important result of Eq. (|77() , whose derivation is one of the aims of this appendix. 

As to the step-like perturbation in this case, we proceed as follows. Let us examine the case (see Eq. HB-7H : 



E{u) 



tp(u). 



1 - ip(u) du 

Using the Taylor series for ip{u), we have (for simplicity we set T = 1) 



(B-29) 



V>N = r(i- M )^- 1 ( At -i) 



,1-/1 



l + u 



r(2 - fi) r(3 - M ) 



(B-30) 



Plugging into Ea. lB-29|l and keeping the lowest orders we obtain: 



E (u) « ^ + — T , (B-31) 



corresponding to 



(B-32) 
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Let us make a check to confirm that E(t) — > /i — 1 for t — > oo. Let us go back from Ea. (|B-29|l to time dominion. 

E(t) -J E(t- t)iP(t)<It = ti(j(t). (B-33) 
o 

So far we did not make any approximation. Let us check if E{t) sa fi — 1 is a solution of Eq. I|B-33(1 . We write 

which in the limit of t — » oo make Eq . (|B-34f) become an identity, thereby proving Eq. (|81|l . 
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(r + iy 



■dr 



(t + l)" ! 



(B-34) 
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